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The diffusion quantum Monte Carlo technique is used to solve the many-body Schrodinger equa- 
tion fully quantum mechanically and nonperturbatively for bosonic atomic gases in cigar-shaped 
confining potentials. By varying the aspect ratio of the confining potential from 1 (spherical trap) 
to 10000 (highly elongated trap), we characterize the transition from the three-dimensional regime 
to the (quasi-)one-dimensional regime. Our results confirm that the bosonic gas undergoes "fermion- 
ization" for large aspect ratios. Importantly, many-body correlations are included explicitly in our 
approach. 



I. INTRODUCTION 



One-dimensional systems of impenetrable bosons are predicted to show intriguing behavior, namely fermionization. 
An impenetrable one-dimensional gas of N bosonic atoms behaves as if the system consisted of ./V one-dimensional 
fictitious spin- less non- interacting fermions, a so-called Tonks or Tonks-Girardeau gas fj}-|l3|. Fermionization has 
profound effects on the characteristics of a system. The ground state energy of N non-interacting bosonic atoms 
in a one-dimensional harmonic trap with trapping frequency v z is E^' z = Nhv z /2. For comparison, the ground 
state energy of N non-interacting spin-less fermionic atoms under the same confinement is Eq' z — N 2 hv z /2, or 
Eq' z = NEq' z . In this paper we map out the transition of a bosonic gas from three- to one-dimensional, and monitor 



the signatures of fermionization |8| jl0| , p.3| 

A Tonks gas has not been observed experimentally yet, and it is an open question how one would probe a Tonks 
gas once it has been created. Effectively one-dimensional atomic gases may be realized by confining N atoms in 
highly elongated cigar-shaped traps where the two tight confining directions are characterized by the frequency v p , 
and the weak confining direction by v z . Correspondingly, the system exhibits two characteristic length scales, the 
oscillator length in the p direction, a p = ykj (muj p ), and the oscillator length in the z direction, a z = y/h/(mu) z ), 
where p = \/ x 2 + y 2 , oj p . z = 2ni' P}Z , and m the atomic mass of the trapped gas. Such a confining potential has been 
realized experimentally by Gorlitz et al. with aspect ratios L = v p jv z as large as 100 @ (see also 0-0). These 
experiments enter the quasi-one-dimensional regime, which is characterized by Hu p 3> fko z , but not the truly one- 
dimensional regime.Quasi-one-dimensional trapping confinement has also been realized utilizing optical lattices fT§| ] 
and wave guides [[l9pc| l. Although experimental techniques will advance further, there will always be an upper limit on 
the ratio a p /a z . The natural question then is under what conditions does the system behave truly one-dimensional, or, 
to which extent is the motion in the tight confinement directions decoupled from the motion in the weak confinement 
direction. 

Various authors have studied the crossover from the three-dimensional regime to the Tonks regime. Common to 
these studies is some sort of variable separation of the tight and loose confinement directions pO| , pl] , p^ 21-24J. For 
a cylindrical trap, e.g., Dunjko et al. [[l0| solve the Schrodinger equation by employing the Lieb-Liniger model j^] 
and a local density approximation. For a toroidal trap, in contrast, Das et al. jK5| develop a variational theory that 
factorizes the wave function into a longitudinal and a transverse part. Our approach here is different in s piri t, in that 
it treats the 3A-dimensional Hamiltonian using the diffusion quantum Monte Carlo (DMC) technique |25|-|2q| fully 
quantum mechanically and nonperturbatively. Thus, no separation of variables needs to be made, and our results 
could be used to benchmark variational calculations. 

The DMC method has been previously applied to microscopic Bose-Einstein condensates (BECs) under spherical 
confinement (29) (for applications of Monte Carlo techniques to gaseous condensates see @-|4| ) • Similarly to Ref. [^9| , 
we monitor the ground state energy and structural properties as a function of the number of particles. In addition, 
we vary the aspect ratio over a wide range, i.e., L — 1, • • • , 10000, which allows signatures of fermionization to be 
investigated. Importantly, our approach includes correlation effects explicitly by treating all 3N coordinates of the 
system equivalently. 

Section |l] outlines the DMC formalism and describes the system under study, while Sec. HI presents our results for 
the energetics and structural properties. Finally, Sec. |^ concludes. 
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II. NUMERICAL TECHNIQUES AND CHARACTERIZATION OF THE SYSTEM 



A. Diffusion quantum Monte Carlo technique 



The DMC method shows a favorable scaling with the number N of particles, and allows bosonic many-body 
systems to be treated fully quantum mechanically and nonperturbatively. The DMC method p5|-p8fl solves the 
time-independent Schrodinger equation 

Hcp n = E n ip n (1) 

for its ground state, n — 0, by considering the time- dependent Schrodinger equation. Introduction of imaginary time 
t = it/% is used to remove the explicit dependence on "i", 

-LT, V ^ r) - [V(f) - E ref ]W, r) . (2) 



dr 2m 



diffusion 



source/sink 



Here, r = {fx,-- - ,tn) collectively denotes the nuclear Cartesian coordinates of the atoms in the trap. V(r) is the 
potential energy surface of the system, which is assumed to be known for now. E re t shifts the absolute energy scale 
without changing the solution to the time-independent Schrodinger equation, Eq. (|l|) . 

A random walk technique is then used to calculate the steady state solution [eigenvector tpo {r) with eigenvalue Eq] 
by propagating to large r: lim,-.^ ^(r, r) — > (fo(r). Typically, about 1000 "walkers" are simulated simultaneously. A 
walker is characterized by its position in Cartesian space and a weight factor indicating its "importance" . To perform 
the propagation in imaginary time r, the short-time approximation to the Green's function is introduced, and the 
kinetic energy part ("diffusion") and the potential energy part ("source/sink") of the Hamiltonian are simulated 
separately at each time step At. 

The DMC algorithm outlined in the previous paragraphs is often refered to as DMC without importance sampling 
or DMC with a constant trial wave function. It turns out that a trial wave function ?At(^) that approximates the 
true ground state wave function t/?o (r) , can be exploited fruitfully for increased accuracy. Introduction of ipT does not 
change the solution to the time-independent Schrodinger equation [Eq. (|IJ)]; however, it does give an additional term 
( "drift" ) in the time-dependent Schrodinger equation, 



dr 2m 



m 



diffusion drift 
-h 2 ( N 



2m 

J =1 



ref 



(*Vt) • (3) 



source/sink 

At each time step Ar, the new drift (or force) term introduces an additional update of the Cartesian coordinates of 
the walkers. Incorporation of a "good" trial wave function tpf, namely one that resembles the true ground state wave 
function (po{r) as closely as possible, leads to a large reduction of the time step error of the DMC calculation [p5|^8[. 
A "bad" trial wave function, in contrast, can lead to biased sampling of the configuration space, and can, in the 
worst case, lead to a seemingly stable solution that differs from the true, stationary solution of the time-independent 
Schrodinger equation, Eq. (nj). Incorporation of a trial wave function ipT is required in the present application to 
obtain converged DMC results. The algorithm with a non-constant trial wave function is typically refered to as DMC 
with importance sampling. 

The DMC algorithm results in an estimate of the ground state energy Eq, which is accompanied by a systematic 
time step error and additional statistical noise due to the stochastic nature of the simulation process. The time step 
error can be reduced by performing a series of calculations for different time steps Ar, and then extrapolating Eq(At) 
to the zero time step limit At = 0. The statistical noise of the energy expectation value can be reduced by performing 
longer simulation runs. 

The calculation of structural ground state properties requires a little more care than the calculation of the ground 
state energy because Eq. (||) interprets the wave function itself as the density of the system, whereas Eq. (||) interprets 
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the product tyipx as the density of the system (see also Sec. [II B| ). To obtain structural properties, such as the radial 
distribution, from the DMC simulation with importance sampling, one therefore either uses approximate mixed 
estimators or employs essentially exact extrapolation techniques such as descendant weighting ]39|-f4l|] . The latter is 
employed in the present application. 



B. Application to trapped atomic gases 



This section discusses the application of the outlined DMC formalism with importance sampling to trapped bosonic 
gases. The DMC method "interprets" the wave function of the system as its density, and its application is thus 
restricted to bosonic systems with a positive and nodeless wave function p2] | . What are the implications for treating 
condensed atomic gases by the DMC method? Ideally, one would like to simulate a condensate using a sum of "exact" 
two-body atom-atom interaction potentials, e.g., the Krauss-Stevens Rb-Rb potential (and possibly three-body 
terms), however, it is well known that the ground state of such a many-body interaction potential corresponds to 
a "metallic snowfiake-like" state. The gaseous condensate state is not the ground state of the system but a highly- 
excited metastable state with a finite lifetime. To treat metastable condensates by the DMC method we artificially 
modify the two-body interaction potential so that the gas-like BEC state becomes the ground state. 

The full potential surface of the N atom system is composed of an internal and an external piece. We parameterize 
the interatomic many-body potential surface by a sum of two-body potentials V2, Ylj>k^( r kj)^ where r^j denotes 
the distance between atom k and atom j, r^j = \rj — r%\. In addition, each atom feels the cigar-shaped external 
confining potential Vt r ap(rk) — 5 TO ( w pPfe + w z z fe)- Our rnany-body Hamiltonian H for TV mass m atoms then becomes 

ff = -2^E V ?+ E M^O+EWr;), (4) 

j—1 j>k,k—l j — 1 

where the Cartesian position vector Tj = (xj, yj, Zj) is measured with respect to the center of the trap. Our calculations 
are performed for a simple two-body hardcore potential , VgjV fej) = 00 for rkj < a, and V2(j*kj) = otherwise, where 
a denotes the two-body s-wave scattering length (2^-||jQ,Q. 

Our trial wave function ipT is chosen to have the simple Bijl-Jastrow-type form ]46[|47|] , 



Mn, ■ ■ ■ ,?*) = ( n ) x ( n *w ) . ( s ) 

which has the proper symmetry for a bosonic system, namely, exchange of two particles leaves ipx unchanged. 4>{ r kj) 
is determined by the solution to the radial two-body scattering equation for V^r^ ). For the hardcore potential, one 
obtains 4>{ r kj) = (1 — o-l^kj) for rkj > a, and 4>(rkj) = otherwise. Thus, <f> prevents two atoms from "crashing 
into each other". In contrast, $ is the part of i/jt that keeps the atoms inside the trap. has two variational 

parameters, p\ and p2, 



$(n) = cxp 



~\ (Pi/ a P f -Pi (zi/a z ) P2 



(6) 



We determine the best parameters p\ andp2 by minimizing the energy expectation value (iPt\H\i}>t) / (iPt^t) through 
variational quantum Monte Carlo (VMC) calculations (2^J|^]. For an ideal bosonic gas (a = 0), p\ — 1/2 and P2 = 2, 
while for a > 0, the optimal values for p\ and P2 depend on a, N, and the strength of the confining potential. 

Our calculations are performed for JV = 2, • • • , 20 with m = m( 87 Rb), v z — 77.78Hz, a two-body s-wave scattering 
length a = 100a. u., and varying v p , v p = \v z , • • • , 10000^ z . The DMC calculations become increasingly demanding 
with increasing aspect ratio L — v p jv z ^ since the strength of the tight confining direction predominently determines the 
timestep At. To ensure an acceptance rate of over 98% in the acceptance/rejection step of the DMC algorithm p7| , 
At has to be decreased with increasing v p as it scales roughly with a p . Thus, to achieve the same total propagation 
time for a system with L = 1 and for a system with L = 10000, we need to perform about two orders of magnitude 
more propagation steps for the system with L = 10000 than for the system with L = 1. 

The number of propagation steps also determines to first order the statistical uncertainty of any quantity calculated 
by the DMC technique. Assume we wish to determine DMC energies with a statistical uncertainty of < 0.2hui z . For 
a bosonic gas with N = 10 atoms, e.g., this implies an accuracy of about 1% for L = 1, and of about 4 x 10 _4 % 
for L — 10000. Achieving convergence of structural properties is more demanding than achieving convergence of 
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the ground state energy. The descendant weighting technique |39|-|41| employed here requires a delay between the 
"original sampling" and the "sampling that is being used for reweighting purposes" . To obtain converged structural 
properties, this delay is as big as 10000 for some of our simulations. 



III. RESULTS 



This section presents our DMC results. Section 



III A 



summarizes energetics, while Sec. IIIBj discusses structural 



properties. Section III C relates our results to those by Petrov et al. M and by Dunjko et al. jl 



A. Energetics 



Figure [T] shows our DMC ground state energies Eq MC (diamonds; for each N, dotted lines connect the DMC 
energies to guide the eye) in units of hcu z as a function of the aspect ratio L on a log-log scale for N = 3, 5, 10 and 20. 
L extends over four orders of magnitude, and Eq MC over about five orders of magnitude. The statistical uncertainty 
of Eq MC is smaller than the symbol size. Figure [j] also shows the ground state energies Eq p (triangles) j|8| obtained 
by solving the mean-field Gross-Pitaevskii (GP) equation [^ in cylindrical coordinates with interaction parameter 
a = (N-l)a, 



11 V72 i i 2 2 I 2 2\ i 



4irh 2 



rr> 



a\xo{p,<P,z)\ 2 



(7) 



and 



The definition of a stems from the derivation of the GP equation based on the Hartree-Fock formalism [p2| 
Xj denote the orbital energy and the orbital wave function, respectively; the GP ground state energy Eq can be 
obtained from the energy functional Eq P [xo] ■ The numerical solution of Eq. (Q) can be simplified by separating out 
the motion in the <p coordinate. On the scale of Fig. [j], the DMC and GP energies are indistinguishable. 



To interpret our results, we subtract the ideal bosonic gas energy in the tight confinement direction, E { 



Nhv n 



from our DMC [GP] ground state energies, and denote the resulting energies by Eq MC,z [Eq P ' z ]. Figure || shows 



E 



DMC,z 




/N (diamonds) as a function of L on a logarithmic scale. To guide the eye, dotted lines connect our data for 
each N. The statistical uncertainty of the DMC energies Eq MC ' z increases with increasing aspect ratio L. E.g., the 
ground state energy of a gas with N — 5 and L = 10 is Eq MC = 52.838(2)ft,i/ z . Subtracting the ideal gas energy in the 
p direction, Eq 
we find Ei 



B,p 



50his z , gives a value of Eq° MC ' z = 2M8(2)hv x with an uncertainty of 0.002/ii/ z . For L = 10000, 
12A(2)hu Z) a 100 times larger statistical uncertainty. Reducing the statistical uncertainty for 
L = 10000 from Q.2hu z to 0.002/ii/ 2 is possible in principle, however, would require significantly more computer time. 

For comparison, Fig. || also shows the GP energies Eq P ' z /N in the z coordinate only (triangles; for each N, dotted 
lines connect these points to guide the eye). For small L, the many-body DMC energies 

e dmc,zj n and the Gp 

energies Eq P ' z /N agree well as expected. However, as L increases, significant devi ations occur. These deviations are 



DMC.z 



a consequence of considering the energy in the z direction only (see also Sec. [II B ); the deviations between the total 
DMC and GP energies, E° MC and £ GP , are negligibly small (see Fig. |l|). As we discuss no w, Fig. ^ confirms that 
the systems under study undergo fermionization for large L as suggested by several authors [^,Eo[^3) . 

Solid horizontal lines on the right hand side of Fig. || indicate the energy in the Tonks limit, E ' z /N = Nhv z /2. For 
each N [[33), our DMC energies Eq M ' z /N approach the Tonks energy, suggesting that the simulation parameters are 
such that the considered systems behave fermionic in the z direction, while the motion in the p direction is frozen out. 
Importantly, our simulations are performed using a trial wave function ipT, Eq. ([b|), with the proper bose symmetry 
for all L. 

To illustrate the fermionization transition further, the lower part of Fig. uiouumo mc ulvlk^ ciragim jj 
by the Fermi energies E P ' Z as a function of L (diamonds) on a logarithmic scale for N = 2, 3, 5, 10 and 20. For each 
N, this ratio is close to zero for small L, and approaches one as L increases. The N — 2 system approaches the Tonks 
limit at smaller L values than the N = 20 system. In addition, the upper part of Fig. || shows the ratio between the 
DMC energies E® M ' z and the ideal bosonic gas energy in the z direction, E^ ' z (triangles). For each N, this ratio is 
close to one at small L, and increases with increasing L. In summary, the energy ratios shown in Fig. || indicate that 
atomic gases under cigar-shaped confinement behave predominantly bosonic for small L and predominantly fermionic 
for large L. 

Our analysis in this section depends profoundly on subtracting the ideal bosonic gas energy in the p direction from 
the many-body energy. The next section discusses structural properties, for which the behavior in the p direction and 



shows the DMC energies E® M ' z divided 
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in the z direction can be naturally separated by integrating over all coordinates but the coordinate of interest, i.e., 
without making any approximations. 



B. Structural properties 



Figure [i] shows the DMC density Uq M ' p , normalized such that / °° n® M ' p (p)dp = 1, for N = 5 for two different 
aspect ratios, i.e., L = 10 (diamonds) and L — 10000 (triangles). The p coordinate is measured in units of a" 1 with 
= 7307a. u. for L = 10, and a p = 231a. u. for L = 10000, respectively (see Tab. |). The statistical uncertainty of 



DMC,p 



is largest for small p since the DMC walk samples the small p region less than the large p region. A solid line 



shows the density for the ideal bosonic gas, 



^(p) = exp [™(p/a p ) 2 ] . (8) 



On the scale of Fig. |J, the DMC densities and the ideal gas density are indistinguishab le, wh ich suggests that the 



motion along the p coordinate is indeed frozen out as L approaches large values (see Sec. Ill A). 

The density along the z coordinate is more interesting than the density along the p coordinate, since its behavior 
changes dramatically as a function of the aspect ratio L. Figure |B| shows the DMC density n^ MC ' z (z), normalized 
such that n$ M ' z (z)dz = 1, for N = 3 for four different aspect ratios L, i.e., L = 10 (diamonds), 100 (crosses), 
1000 (triangles), and 10000 (stars), z is measured in units of aj 1 , where a z = 23107a. u.. For each L, dotted lines 
connect the symbols to guide the eye. Additionally, a solid grey line shows the one-dimensional ideal bosonic gas 
density n^' z , 



B.z 



{z) = exp [-{z/a z f} . (9) 



A solid dark line shows the ideal fermionic gas density rig' z , which is given by the sum of squares of the single-particle 
wave functions, 

1 1 N ^ 1 

= nT^oZ £ 2^ 2(z/a * )exp H-K) 2 ] . (io) 

k— 

where the denote Hermite polynomials. The Fermi density shows three distinct maxima, corresponding to N — 3. 

To ensure an unbiased observation of the fermionization effect in our DMC calculations, we choose a "structureless" 
trial wave function ipj*, that is a trial wave function without oscillations. Figure |^ shows the z component of "0T 
squared (grey solid line), 

Mexp[-2 Pl (z/a z ) P2 ], (11) 

with pi = 0.018 and P2 = 4.43 (M is a properly choosen normalization constant), together with the DMC density for 
L = 10000 and a = 100a. u. (diamonds) and the ideal fermionic gas density [Eq. (|l0|), black solid line]. We checked 



the "delay" employed in the descendant weighting scheme of our DMC calculation (see Sec. II B), and estimate the 
statistical uncertainty to be about three times the symbol size. The DMC density for L = 10000 does not seem to 
show any Fermi-oscillations (note, however, that the statistical error of the DMC density is non-negligible, as indicated 
exemplary for one data point at z = la^T 1 ). 

Based on a single-particle picture, Petrov et al. || and Dunjko et al. |]Io| derive an expression for the density of a 
gas under cylindrical confinement in the Tonks limit, 



^ = ^^|(, 1 -2^f) 1/2 - ^ 

Notably, nJ j onks is dependent on N and a z only, i.e., it is independent of a. Figure^ indicates that n^ onks (dashed 
line) follows the Fermi density [Eq. (|l0|), dark solid line] closely, except for the oscillations, which cannot be described 



within a single-particle picture. The DMC density n ' agrees well with . Section [II C presents a detailed 

comparison between our work and that by Petrov et al. §| and by Dunjko et al. |]l0|. 
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The Fermi density, Eq. (|To|), changes significantly with increasing N. To illustrate this effect, Fig. ^ shows the 
DMC density n^ MC ' z as a function of z for N = 10 for three different L values, L = 10 (diamonds), 100 (crosses), 
and 1000 (diamonds), using the same scale as in Fig. ^| for N — 3. The Fermi density (dark solid line) now shows 
ten maxima, owing to the increased number of terms in Eq. ([l0]) . As L increases the DMC density approaches the 
Fermi density. Our density s tudie s, discussed here in detail for N — 3 and N — 10, support our interpretation of the 
energetics presented in Sec. [II A, and show strong evidence that an atomic gas under external confinement indeed 
undergoes fermionization as L increases. 



C. Comparisons with the literature 



This section relates our results to a set of criteria put forward by Petrov et al. M and by Dunjko et al. [JlO|. Dunjko 
et al. |Tc(l derive five requirements for a successful experimental observation of the zero-temperature Tonks gas density 
profile, Eq. (HI). 



1. a positive and CI — < 1, where C as 1.4603. These conditions imply a positive one-dimensional coupling 

V 2 p 

constant gid, 



gid 



2Ti 

mai d 



with aid 



a K l _ 



V2% 



(13) 



2. C2 



<C 1. The transverse energy level spacing must exceed the energy of the gas in the z direction. 



3. C3 = T) -C 1, where rj denotes the governing parameter of the system, 



7] = 



(9/2) 



1/3 



4 a 4 I _ C_JL 



2/3 



N 2/3 . 



(14) 



4. C4 = 2ti i ?(ma i ) ^ ener §y ma' 1 must exceed the energy level spacing in the z direction. 

' > Id' Id 

5. N > 1. This condition is fulfilled trivially. 

Table || summarizes the characteristic lengths a p and the one-dimensional scattering lengths aid as a function of L 
for a — 100a. u., m = m( 87 Rb), and v z — 77.78Hz. Condition 1. puts an upper limit on the aspect ratio L; for the 
parameters chosen here, the coupling constant gid becomes negative (aid becomes positive) for L a 50000. CI and 
C4 are independent of N, while C2 and C3 scale with N and iV 2 / 3 , respectively. 

Figure || shows the criteria CI through C4 as a function of L on a log-log scale. C2 and C3 are shown for N = 3, 
and 20. C2 — C4 are smaller than 0.1 for about L > 3500, and smaller than 0.01 for the extreme aspect ratio of about 
L > 20000. The behavior of C2 and C3 indicates that the "window" , in which an atomic gas behaves like a Tonks gas, 
becomes narrower as N increases. It thus seems advantageous to experimentally work with small particle samples. 
Toward this end, optical lattices |l^,Q, which allow preparation of several 1000 identical copies of microscopic atom 
samples, seem most promising. Additionally tuning the two-body s-wave scattering length to large positive values 
through the use of Feshbach resonances |^5|,[56| would move the transition to the Tonks regime to smaller aspect 
ratios. A critical question, which is beyond the scope of this paper, is whether the density profile in the z direction 
can be unambiguously measured in time- of- flight expansion experiments, or whether mean- field effects may mask the 
signature of the Tonks density profile p2| , p3|j57| ] . 



IV. SUMMARY AND OUTLOOK 



The present DMC studies of atomic gases under external confinement with varying aspect ratio provide first bench- 
mark results, which will allow the accuracy of numerically less demanding, and more commonly used approximate 
approaches to be tested. In contrast to variable separation schemes and mean-field-type approaches, the DMC for- 
malism treats many-body correlations explicitly. The DMC method consistently treats atomic gases under spherical 
confinement and under highly-elongated cylindrical confinement, and thus allows the transition from three-dimensional 



G 



to one-dimensional to be mapped out uniquely. Application of the DMC technique can be extended to arbitrary con- 
fining potentials such as toroidal or random potentials. The main modification of the DMC algorithm would amount 
to adjusting the $ component of the trial wave function, Eqs. (JsJ) and accordingly. 

The DMC method in its present implementation is restricted to treating many-body systems in their ground state. 
To nevertheless treat metastable atomic gases, the present study parameterizes the atom-atom interaction potential 
by a simple repulsive model potential, i.e., a hardcore potential. It would be interesting to investigate dependencies 
of the energetics and structural properties on the detailed shape of the two-body interaction potential as a function 
of the aspect ratio L. Of particular interest would be the characterization of atomic gases under external confinement 
as a function of the ratios a/a z and b/a z , where b denotes the characteristic length of the two-body potential [^8[-|60[. 

The application of the zero-temperature DMC technique in its present implementation is restricted to a description 
of static properties: We cannot make predictions regarding the lifetime or stability of trapped atomic gases. Our studies 
here characterize atomic gases, assuming that they are stable. Our approach provides no access to information related 
to finite temperature effects, time evolution, and coherence/decoherence phenomena. 

In summary, this paper presents exact many-body calculations for trapped atomic bosonic gases under external 
confinement with varying aspect ratio. The trapping geometries considered range from spherical to highly elongated, 
L = 1 — 10000. The DMC technique is powerful in assessing the importance of many-body correlations. Our study 
provides first benchmark results for atomic gases under highly elongated cylindrical confinement. 
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FIG. 1. Ground state DMC energy, E° MC (diamonds; 
for each TV, dotted lines connect the DMC energies to guide 
the eye), together with the ground state energy obtained by 
solving the two-dimensional Gross-Pitaevskii equation, Eq P 
(triangles), as a function of the aspect ratio L = v p /v z on a 
log-log scale for N = 3,5, 10 and 20. The statistical errorbar 
of the DMC energies is smaller than the symbol size. 



10 














DMC 

GP 

Ideal 






N=20 . 


8 




Id FG 




6 








10 i 

1 b 




4 






N=20 a 

to 


/ ff 


N=10 . 


2 






A/ 


A ^ .© <>♦««»> 


N=5 _ 







1 




N=3 



10° 10 1 10 2 10 3 W A 
aspect ratio L=u p /u z 



10" 



FIG. 2. Energy per particle along the z coordinate only. 
The DMC energies Eq MO ' z /N (diamonds) are shown to- 
gether with the GP energies E° p ' z /N (t riangles) as a function 
of the aspect ratio L on a logarithmic scale for TV = 3, 5, 10 
and 20. To guide the eye, dotted lines connect the DMC 
and GP energies for each TV. Clear discrepancies between 
Eq MC ' z /N and Eq P ' z /N are visible at larger L. Horizontal 
solid lines on the right hand side show the energies in the 
Tonks limit, E F ' Z /N. 
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FIG. 3. Ratio between the DMC energy in the z direction 
and the Tonks energy, Eq MC ' z / E F ' Z (lower part, diamonds) 
as a function of L on a logarithmic scale, together with the 
ratio E° mc ' z /Eq' z (upper part, triangles) for TV = 2,3,5, 10 
and 20. As L increases the ratio Eq MC,z /E f ' z approaches 
one, while the ratio Eq MC ' z /Eq' z approaches values much 
greater than one, i.e., the value TV, indicating that the con- 
fined atomic gas undergoes fermionization as L increases. 
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4. DMC density n ' p as a function of p for 
for two different aspect ratios L, L = 10 (diamonds) 
= 10000 (triangles). To collapse the two DMC densi- 



ties onto the same scale, p is measured in units of a p with 
a p = 7307a.u. for L = 10, and a p = 231a.u. for L = 10000, 
respectively. Additionally, a solid line shows the ideal bosonic 
gas density [Eq. ||]. The DMC densities follow the ideal 
bosonic gas density closely. 
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FIG. 5. DMC density n° MC ' z for TV = 3 as a function 
of z for four different aspect ratios, L — 10 (diamonds), 100 
(crosses), 1000 (triangles), and 10000 (stars), z is measured 
in units of a" 1 , where a z — 23107a.u.. In addition, a grey 
solid line shows the ideal gas density for bosons [Eq. (HI)], and 
a dark solid line shows that for fermions [Eq. (|lo|)] . The DMC 
density is close to n^' z for small L, and approaches n F ' z as L 
increases. 
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FIG. 6. DMC density °' z for AT = 3 as a function 
of z for L = 10000 (diamonds; the statistical uncertainty is 
indicated exemplary for z = laj 1 ). z is measured in units 
of a' 1 with a z = 23107 a.u.. In addition, a dark solid line 
shows the fermionic density n^' z (z) [Eq. (|lo|)], and a grey 
solid line the square of the trial wave function component in 
the 2 direction [Eq. (|TT[)] . The dashed line shows the density 
profile for the Tonks limit [Eq. (jT|)]. 
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FIG. 8. Criteria CI (grey solid line), C2 (grey dotted 
hnes), C3 (dark dashed lines), and C4 (dark solid line) as 
a function of the aspect ratio L on a log-log scale usin g 
a = 100a.u., m = m( 87 Rb) and v z = 77.78Hz (see Sec. |lIIC| ). 
CI and C4 are independent of N, while C2 and C3 scale with 
N and N 2 ^ 3 , respectively. C2 and C3 are shown for N — 3 
and N = 20. 
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FIG. 7. DMC density n° MC ' z for JV = 10 as a function 
of z for four different aspect ratios, L — 10 (diamonds), 100 
(crosses), and 1000 (triangles), z is measured in units of a^ 1 , 
where a z — 23107a.u.. In addition, a grey solid line shows the 
ideal gas density for bosons [Eq . (S)], and a dark solid line 
shows that for fermions [Eq. The DMC density is close 

to Uq ,z for small L, and approaches n^' z as L increases. 



L 


a p [a.u.] 


aid [a.u.] 


1. 


23107. 


-5.3 x 10 b 


5. 


10328. 


-1.1 x 10 6 


10. 


7307. 


-5.3 x 10 5 


50. 


3268. 


-1.0 x 10 5 


100. 


2311. 


-5.1 x 10 4 


500. 


1033. 


-9.6 x 10 3 


1000. 


731. 


-4.6 x 10 3 


5000. 


327. 


-7.3 x 10 2 


10000. 


231. 


-3.0 x 10 2 


50000. 


103. 


-8.3 x 10 -2 



TABLE I. Characteristic length a p and one- dimensional 
scattering length aid f° r a = 100a.u., m = m( 87 Rb), 
v z — 77.78Hz, and a z = 23107a.u. as a function of the aspect 
ratio L. 
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